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ABSTRACT 



After new observations of 39 galaxies at z ~ 0.6-1.0 obtained at the IRAM 30m telescope, we present our full CO line survey covering 
the redshift range 0.2 < z < 1. Our aim is to determine the driving factors accounting for the steep decline of the star formation rate 
during this epoch. We study both the gas fraction, defined as M gal /(M gal +M IM ,-), and the star formation efficiency (SFE) defined by the 
ratio between far-infrared luminosity and molecular gas mass (Lpi R /M(H 2 )), i.e. a measure for the inverse of the gas depletion time. 
The sources are selected to be ultra-luminous infrared galaxies (ULIRGs), with Lfir larger than 10 12 L G , experiencing starbursts. 
Adopting a standard ULIRG CO-to-H 2 conversion factor, their molecular gas depletion time is less than 100 Myr. Our full survey 
has now filled the gap of CO observations in the 0.2<z<l range covering almost half of the cosmic history. The detection rate in 
the 0.6 < z < 1 interval is 38% (15 galaxies out of 39), compared to 60% for the 0.2<z<0.6 interval. The average CO luminosity 
is L'co = 1.8xl0 10 K kms~' pc 2 , corresponding to an average H2 mass of 1.45x10'° M Q . From the observation of 7 galaxies in 
both CO(2-l) and CO(4-3), a high gas excitation has been derived; together with the dust mass estimation, this supports the choice 
of our low ULIRG conversion factor between CO luminosity and H2 for our sample sources. We find that the gas fraction in our 
sample increases strongly with redshift, and is the most determining factor in the cosmic star formation evolution. The SFE evolution 
in comparison is shallower, and we conclude that it plays a less significant role in setting the cosmic star formation evolution. 

Key words. Galaxies: high redshift — Galaxies: ISM — Galaxies: starburst — Radio lines: Galaxies 



1. Introduction 

The star formation history (SFH) of the Universe shows a steep 
decline by a factor 10 between z= 1 and 0, after a peak of activity 
around z=l-1.5, well within the first half of the Universe's age 
(Madau et al 1998, Hopkins & Beacom 2006). During the sec- 
ond half of the Universe, fundamental changes occur not only in 
the star formation rate, but also in the galaxy interaction/merger 
rate (Le Fevre etal. 2000, Conselice etal. 2009, Kartaltepe 
etal. 2010), and in the galaxy morphology (Sheth etal. 2008, 
Conselice et al. 201 1), as the well known local Hubble sequence 
is just assembling. 

The main physical processes causing the decline in the SFH 
are not well understood. One frequently invoked factor is the gas 
content, or gas fraction at a given stellar mass, since gas is the 
fuel for star formation. Locally, the gas fraction for giant spirals 
is about 7-10% (Leroy et al. 2008, Saintonge etal. 201 la), while 
at z~1.2 it increases to 34+5% and at z~2.3 to 44+6% (Tacconi 
etal. 2010, Daddi etal. 2010). However, the average surface and 
volume gas densities should be more relevant factors, and there 
is evidence of a wide range in gas properties of galaxies at each 
redshift: molecular gas extents, mean gas densities and star for- 
mation efficiencies (Daddi etal. 2008). The average molecular 
gas density can be traced by the CO excitation, measured from 
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line ratios of the CO ladder (e.g. Weiss etal. 2007). Another in- 
tervening factor is an external dynamical trigger of star forma- 
tion, such as galaxy interactions, or accretion by cold gas. Both 
effects are expected to increase with redshift, up to z=2-3. 

Galaxy mergers provide the violent gravity torques able 
to drive the gas quickly to galaxy centers, and trigger star- 
bursts. Locally Ultra-Luminous Infra-Red Galaxies (ULIRGs), 
with Lfir > 10 12 L Q , are in great majority starbursts caused by 
galaxy major mergers (e.g. Sanders & Mirabel 1996, Veilleux 
etal. 2009). The global gas fraction in those systems is not high, 
and they are expected to sustain such a high rate of star for- 
mation only for a short duration, on the order of lOOMyr. The 
depletion time-scale of massive star-forming galaxies, with high 
gas fraction but low star formation efficiency (SFE), is on the 
other hand of the order of one Gyr, or more (Leroy et al. 2008, 
Bigiel et al. 2008). These galaxies are dubbed to fall on the "main 
sequence" (as defined below) of star forming galaxies (Noeske 
et al. 2007), or in a "normal" star formation phase. 

These two modes of star formation: main sequence (MS) and 
merger-induced starbursts (SB) have been been discussed as a 
function of stellar mass. At both high- and low-z several studies 
(Noeske 2009; Rodighiero etal. 2011; Wuyts etal. 2011 ) have 
shown a broad correlation between the star formation rate (SFR) 
and the stellar mass (M»). The bulk of star forming galaxies fol- 
low the "normal" mode of star formation in the center of corre- 
lation, and define the SFR-M* "main sequence". The galaxies in 
the upper envelope of this main sequence have higher SFR at a 
given mass. These starbursts are rather rare, and represent only 
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10% of the cosmic SFR density at z ~ 2 (Rodighiero et al. 201 1). 
They could be either a common phase of "normal" galaxy life, 
or the end phase leading to a massive early-type galaxy, where 
star formation is quenched. 

The specific star formation rate (i.e. star formation per unit 
mass) sSFR, is found either to decrease slightly with M* (oc 
AC - 2 , Noeske etal. 2007), or to be a constant with M» (Wuyts 
etal. 201 1). The global SFR decreases with cosmic time as SFR 
oc (1+z) 2 - 7 , out to z=l-2. At each epoch, there exists a popula- 
tion of non star-forming (quiescent) galaxies, of high mass, and 
high Sersic index. These galaxies might correspond to quenched 
galaxies and their average stellar mass increases with z. The 
starbursts at the top of the main sequence have a morphol- 
ogy intermediate between the MS and quiecent objects (Wuyts 
etal. 2011). 

To better clarify the role of each physical process in the SFH 
decline in the second half of the Universe's history, it is of prime 
importance to determine the gas content, spatial extent and sur- 
face density of star forming objects between redshifts 0.2 and 

1 . To do so we have begun a CO line survey of ULIRGs in this 
still unexplored redshift range (e.g. Combes etal. 2011, here- 
after paper II, where we focused on objects in the redshift range 
0.2<z<0.6), and we now complement it with 0.6 < z < 1 objects. 

Very little is known about the molecular gas content of galax- 
ies in this redshift range because of observational difficulties. 
Paradoxically, it is often easier to study higher redshift galax- 
ies, as the global CO line flux increases almost like the square 
of frequency for high-J lines (Combes etal. 1999). In the most 
favorable 3mm atmospheric window, between 81 and 115 GHz, 
all redshifts can be observed with at least one line of the CO ro- 
tational ladder, except between z=0.4 and 1 . For galaxies in this 
range, CO lines must be searched for in the 2mm or 1mm range, 
and several CO lines are out of reach. 

In this paper we focus on a sample of 39 IR-luminous galax- 
ies in the range 0.6 < z < 1, detected either by IRAS or Spitzer. 
From their stellar mass and derived sSFR, these objects are 
mostly in the starburst phase on the upper envelope of the main 
sequence. Their sSFR helps us to determine which CO-to-H2 
conversion ratio to use, and therefore provides the best esti- 
mates of the true star formation efficiency. There appears to 
be a variation in conversion ratios by a factor ~ 6 between the 
two extreme classes of objects. In local ULIRGs where a com- 
pact nuclear starburst is typically mapped, the CO emission is 
much stronger for a given H2 mass, and the generally adopted 
ratio between M(H2) and L'co, expressed in units of M (K 
kms" 1 pc 2 )~', is ff=0.8 instead of the adopted Milky Way ra- 
tio of ff=4.6 (Solomon etal. 1997). However, it is still possible 
that this ratio underestimates the molecular mass in ULIRG (e.g. 
Papadopoulos et al. 2012). 

To determine the extent of the molecular gas, interferometric 
maps are required, and we are also mapping some of the de- 
tected galaxies with the IRAM Plateau de Bure interferometer 
(PdBI) (Combes et al. 2006, paper I). This first map of a galaxy 
at z=0.223 revealed that at least half of the emission is extended 
on scales of 25-30kpc. We will present other PdBI maps in a fu- 
ture article, the focus of the present paper is the second part of 
our survey. The sample is described in Sect. [2] and the observa- 
tions in Sect. [3] Results are presented in Sect.|4]and discussed in 
Sect. 

2. The sample 

We used the same criteria as in paper II (except for the redshift 
range), and selected ULIRG galaxies with log Lfir/L q > 12.45, 
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Fig. 1. Redshift histogram of the presented sample (see Table 
Q} at redshifts between 0.6 < z < 1.0. The sources have been 
selected based on their 60 /im flux from IRAS, or 70 //m 
flux from Spitzer-MIPS. The open histogram indicates CO- 
undetected sources, while the filled histogram shows detections. 



between 0.6 < z < 1, with a declination larger than -12°, a spec- 
troscopic redshift, and detected at 60/vm (IRAS), or at 70yimi 
(Spitzer). Our initial sample during the first run contained 39 ob- 
jects; however, observations of some galaxies could not be used 
because of baseline oscillations caused by their strong contin- 
uum flux, or due to noisy spectral parts of the receivers. We have 
therefore removed those galaxies from our sample and replaced 
them with slightly lower far-infrared luminosities, but keeping 
the ULIRG (Lfir > 10 12 L G ) criterion. Among the 39 sources, 32 
have SDSS images; they often appear to be point sources, and in 
a few cases show perturbed and interacting systems. Ten sources 
have been imaged in the near-infrared with sub-arcsec seeing 
by Stanford etal. (2000). HST images are available for about 
half of the sample (17 objects). From the available imaging, it is 
possible to identify morphological perturbations or minor com- 
panions in about 20% of the cases, and obvious strong interac- 
tions in about 30%. The remaining half of the objects appears 
unperturbed. This high frequency of unperturbed morphologies, 
in comparison with our sample between 0.2 < z < 0.6 is at least 
in part due to the lack of high spatial resolution and sensitivity 
optical imaging for the higher-z sample. 

The name of the sources, their coordinates and redshifts are 
given in Table Q] Out of the 39 objects in the sample, 15 were 
detected, corresponding to a detection rate of 38%. Figure[T]dis- 
plays their distribution in redshift. 

The majority of the sources have IRAS fluxes at 60 and 
lOOjurn, and their far-infrared fluxes Ffir are computed as 1.26 
x 10~ 14 (2.58 Sso+Sioo) W nr 2 (Sanders & Mirabel 1996). The 
far-infrared luminosity between 40 and 500jt/m is then Lpm = 
4 n D 2 L CC Ffir, where is the luminosity distance, and CC 
the color correction, CC=1.42. Some of the sources have instead 
Spitzer-MIPS fluxes at 70 and 160yum. We chose to compute 
their Ffir fluxes similarly, with a linear combination of these 
two fluxes, since all sources are ultra-luminous objects, i.e. their 
SEDs fall into the same category and this gives results compara- 
ble to those of SED fitting (e.g. Symeonidis etal. 2008). The 
FIR-to-radio ratio ? =log([F F iR/(3.75 10 12 Hz)]/[f v (1.4 GHz)]) 
has been computed for sources where radio data are available; 
the radio fluxes are listed in Table [2] Excluding the radio galax- 
ies 3C280 and 3C289 (SI 8 and S20), the average is 0=1.7, close 
to the nominal value for ULIRGs (Sanders & Mirabel 1996). The 
star formation rates are all above 200 M yr~' , and the sample av- 
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Velocity (kra/s) Velocity (km/s) 

Fig. 2. CO spectra of the detected galaxies. Zero velocity cor- 
responds to the optically determined redshift, listed in Table Q] 
According to the redshift, the lowest frequency observable is ei- 
ther the CO(2-l) or CO(3-2) line. One source (S16) is detected 
in CO(4-3) but not in CO(2-l), as indicated in Tabled The ver- 
tical scale is T m t, in mK. 



erage is 1200 M yr , as estimated from the infrared luminosity 
(e.g. Kennicutt 1998). 

In this article, we adopt a standard flat cosmological model, 
with A = 0.73, and a Hubble constant of 71 kms -1 Mpc~' 
(Hinshaw et al. 2009). 



3. Observations 

The observations were carried out with the IRAM 30m tele- 
scope at Pico Veleta, Spain, in two periods, January and June 
2011, with a few remaining sources observed in the pool be- 
tween November 2011 and February 2012. According to their 
redshifts, sources were observed in their CO(2-l), CO(3-2) or 
CO(4-3) lines, either with the 3mm, 2mm or 1mm receivers. 
Often only one line was observable. When possible and when 
the atmosphere was favorable, we simultaneously observed two 
CO lines. 

The broadband EMIR receivers were tuned in single side- 
band mode, with total bandwidth of 4 GHz per polarisation. 
This covers a velocity range of ~ 12,000 kms" 1 at 3mm, and 
~ 8,000 kms -1 at 2mm. The observations were carried out in 
wobbler switching mode, with reference positions offset by 2' 
in azimuth. Several backends were used in parallel, the WILMA 
autocorrelator, with 2 MHz channel width, covering 4x4 GHz, 
and the 4 MHz filterbanks, covering 2x4 GHz. 




-2000 -1000 1000 80002000 -1000 1000 2000 
Velocity (km/s) Velocity (km/s) 



Fig. 3. Same as Fig.|2]for additional galaxies. 



We spent on average 3 hours on each galaxy, and reached 
a noise level between 0.7 and 2 mK (main beam temperature), 
smoothed over 30 km s _1 channels for all sources. Pointing mea- 
surements were carried out every 2 hours on continuum sources 
and the derived pointing accuracy was 3" rms. The tempera- 
ture scale used is main beam temperature r m b- At 3mm, 2mm 
and 1mm, the telescope half-power beam width is 27", 17" 
and 10", respectively. The main-beam efficiencies are rj m \, = 
T*J r mb =0.85, 0.79 and 0.67, respectively, and S/T mb = 5.0 Jy/K 
for all bands. Spectra were reduced with the CLASS/GILDAS 
software, and the spectra are smoothed up to ~ 60 km s _1 chan- 
nels for the plots. 

4. Results 

4.1. CO detection in z=0.6-1.0 ULIRGs 

Figures [2] [3] and [4] display the CO-detected sources, in their 
lower-J CO spectrum. For three galaxies, also their higher-J lines 
are plotted. Table [2] reports all line parameters, and also the 
upper-limits for the non-detections. Integrated signals and ve- 
locity widths are computed from Gaussian fits. These give also 
the central velocities, with respect to the optical redshift of Table 
Q] There are some offsets between the CO and optical velocities, 
always lower than 700 km s _1 . These offsets might be due to the 
accuracy of the optical determination, or be intrinsic due to an 
outflow of ionised gas. The upper limits are computed at 3cr, 
assuming a common line-width of 300kms _1 , and getting the 
rms of the signal over 300kms _1 . Lines are considered detected 
when the integrated signal is larger than 3<x. All spectra above 
3cr are shown in the Figures. 
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Table 2. Observed line parameters 



Gal Line v obs S(CO)" V b AVWhm L' co /10 10 S 60 S^ 00 log L F1R F(1.4GHz) rf log M{ 

[GHz] [Jykms- 1 ] [kms- 1 ] [kms- 1 ] [K kms" 1 pc 2 ] [Jy] [Jy] [L ] [mJy] [M Q ] 



SI 


CO(21) 


133.490 


2.1 ±0.3 


354. 


± 37. 


469. 


±73. 


1.53 


0.26 


0.44 


13.09 


4.3 


10.84 


SI 


CO(43) 


266.960 


4.1 ± 1.1 


428. 


± 96. 


702. 


± 208. 


0.74 


0.26 


0.44 


13.09 


4.3 


10.84 


S2 


CO(21) 


143.014 


< 1.2 










< 0.6 


0.23 


0.21 


12.77 


9.0 


10.58 


S3 


CO(21) 


140.145 


< 1.0 










< 0.6 


0.14 


0.03 


12.51 


70.3 


9.58 


S4 


CO(43) 


284.242 


< 1.0 










< 0.1 


0.06 


<0.12 


12.32 


— 


11.55 


S5 


CO(43) 


256.419 


< 3.0 










< 0.7 


0.21 


0.40 


13.12 


2.0 


10.84 


S8 


CO(21) 


139.196 


< 1.2 










< 0.7 


0.03 


0.36 


12.58 


73.7 


10.71 


S9 


CO(21) 


141.957 


< 1.0 










< 0.5 


0.25 


0.03 


12.72 


2.6 


10.20 


S10 


CO(21) 


141.698 


1.3 ± 0.2 


132. 


± 27. 


291. 


± 64. 


0.68 


0.15 


0.03 


12.51 


1.1 


10.59 


Sll 


CO(43) 


250.566 


< 4.1 










< 1.0 


0.25 


0.03 


13.03 


9.3 


11.01 


S12 


CO(32) 


190.837 


4.8 ± 1.4 


20. : 


t 36. 


271. 


± 102. 


1.93 


0.15 


0.03 


12.79 


- 


10.42 


S13 


CO(21) 


142.220 


< 0.8 










< 0.4 


0.22 


0.32 


12.83 


0.7 


11.43 


S14 


CO(21) 


141.175 


1.1 ± 0.2 


1. 4 


: 15. 


122. 


± 29. 


0.61 


0.22 


0.54 


12.94 


6.3 


10.62 


S14 


CO(43) 


282.327 


3.7 ± 0.9 


-82. 


+ 65. 


493. 


± 99. 


0.49 


0.22 


0.54 


12.94 


6.3 


10.62 


S16 


CO(21) 


116.704 


< 3.2 










< 4.2 


0.20 


0.41 


13.33 


4.5 


12.66 


S16 


CO(43) 


233.389 


10.3 ± 1.6 


521. 


± 31. 


483. 


± 103. 


3.36 


0.20 


0.41 


13.33 


4.5 


12.66 


S18 


CO(43) 


230.982 


< 2.0 










< 0.7 


0.12 


0.06 


12.95 


5100.0 


10.33 


S20 


CO(43) 


234.340 


< 1.8 










< 0.6 


0.10 


0.03 


12.82 


2400.0 


11.66 


S21 


CO(21) 


133.614 


<0.8 










<0.6 


0.08 


0.15 


12.60 


5.9 


11.54 


S22 


CO(21) 


138.047 


1.8 + 0.4 


-656. 


± 25. 


249. 


± 57. 


1.11 


0.19 


0.03 


12.68 


11.0 


10.04 


S22 


CO(43) 


276.072 


6.0 + 0.9 


-685. 


± 16. 


214. 


± 33. 


0.91 


0.19 


0.03 


12.68 


11.0 


10.04 


S23 


CO(21) 


140.572 


< 1.0 










< 0.6 


0.28 


0.03 


12.79 


- 


10.29 


S25 


CO(43) 


239.415 


< 3.0 










<0.9 


0.28 


0.50 


13.40 


5.9 


11.36 


S26 


CO(21) 


137.882 


< 1.2 










< 0.7 


0.15 


0.03 


12.59 


1.9 


11.85 


S27 


CO(32) 


181.330 


1.5 ± 0.3 


-128. 


± 28. 


211. 


± 60. 


0.77 


0.03 


0.26 


12.82 


0.2 


11.25 


S28 


CO(43) 


241.383 


< 2.2 










< 0.6 


0.03 


0.24 


12.80 


- 


10.14 


S29 


CO(21) 


134.818 


4.4 ± 0.5 


68. : 


t 18. 


373. 


±48. 


3.02 


0.17 


0.21 


12.84 


1.2 


11.39 


S29 


CO(43) 


269.614 


8.4 ± 1.6 


40. : 


t 38. 


421. 


± 85. 


1.43 


0.17 


0.21 


12.84 


1.2 


11.39 


S30 


CO(21) 


136.494 


2.4 ± 0.4 


-133. 


± 21. 


267. 


± 43. 


1.52 


0.24 


0.23 


12.92 


2.3 


11.04 


S30 


CO(43) 


272.967 


1.3 ± 0.4 


-32. 


± 6. 


44. 


± 26. 


0.21 


0.24 


0.23 


12.92 


2.3 


11.04 


S32 


CO(43) 


240.251 


< 1.8 










<0.5 


0.40 


0.47 


13.48 


- 


10.72 


S34 


CO(21) 


133.259 


3.6 ±0.5 


-60. 


+ 41. 


569. 


± 96. 


2.58 


0.19 


0.56 


13.07 


2.7 


10.84 


S34 


CO(43) 


266.498 


4.8 ± 1.1 


-34. 


+ 52. 


414. 


± 104. 


0.87 


0.19 


0.56 


13.07 


2.7 


10.84 


S35 


CO(21) 


142.022 


1.8 ±0.5 


93. : 


t 19. 


182. 


± 70. 


0.93 


0.07 


0.16 


12.41 


166.4 


11.79 


S36 


CO(43) 


238.881 


< 3.6 










< 1.1 


0.03 c 


0.14 e 


12.66 




11.62 


S37 


CO(43) 


249.211 


< 2.8 










< 0.7 


0.01 c 


0.06 c 


12.21 


~_ 


11.29 


S38 


CO(21) 


138.878 


< 1.4 










<0.8 


0.02" 


0.08" 


12.04 




10.83 


S39 


CO(21) 


130.247 


< 1.2 










< 1.0 


0.01" 


0.10" 


12.25 




12.19 


S40 


CO(43) 


235.225 


< 1.8 










<0.6 


0.01" 


0.09" 


12.42 




11.69 


S41 


CO(21) 


136.413 


< 1.0 










<0.7 


0.01" 


0.1 2 e 


12.15 




11.77 


S42 


CO(21) 


143.472 


2.3 ± 0.4 


-98. 


+ 93. 


1001. 


± 244. 


1.12 


0.06 


0.20 


12.41 


15.0 


11.86 


S43 


CO(21) 


138.878 


4.7 ± 0.6 


560. 


±23. 


385. 


±49. 


2.75 


0.05" 


0.14 e 


12.38 




11.05 


S44 


CO(21) 


138.047 


5.9 ± 1.0 


60. : 


t36. 


415. 


± 72. 


3.59 


0.07" 


0.25" 


12.57 




11.38 


S45 


CO(21) 


132.345 


< 1.8 










< 1.4 


0.04" 


<0.15 e 


12.47 


1.2 


11.67 


S46 


CO(21) 


138.212 


< 1.4 










<0.9 


0.04" 


<0.16 e 


12.39 


4.6 


10.72 


S47 


CO(21) 


131.030 


3.1 ±0.7 


46. : 


t46. 


437. 


± 153. 


2.45 


0.03 e 


<0.11 c 


12.34 


0.5 


11.31 



Quoted errors are statistical errors from Gaussian fits. The systematic calibration uncertainty is 10%. 
" The upper limits are at 3cr with an assumed AV = 300 km s" 1 , except for S 16, where AV is known. 
* The velocity is relative to the optical redshift given in Table[T] 

The 60 and 100/jm fluxes are from NED (http://nedwww.ipac.caltech.edu/ 1 or Stanford et al. (2000) 

d From the FIRST catalog (http://sundog.stsci.edu/ 1. Errors are typically 0.14 mJy 

' For these objects, the Spitzer 70 and 160/jm fluxes replace the IRAS 60 and 100/iin ones. 

' Stellar masses were obtained through SED fitting (see section l4~6t . 



The detection rate of 38% in this z=0.6-l sample is lower 
than the value of 60% found in the lower-z counterpart at z=0.2- 
0.6 (paper II), however this can already be explained by the re- 
duced sensitivity for the more distant objects. Indeed, the lines 
detected are now CO(2-l) instead of CO(1-0), and for excited 
point sources the flux could be up to ~ 4 times larger, however, 
the average luminosity distance at 0.6 < z < 1.0 weakens the 



signal by about a factor 5.5 with respect to the previous closer 
sample, and this outweighs the gain by climbing up the CO lad- 
der. At 0.6 < z < 1.0, the average angular distance is 1550 Mpc, 
i.e. the beam of 17" corresponds to 127 kpc, so the sources can 
all be considered as point-like, at this resolution. 

The line-widths detected are compatible with what is ex- 
pected from massive galaxies of ULIRG type, and quite com- 
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Table 1. Definition of the sample, selected by FIR luminosity 




-1000 1000 

Velocity (kra/s) 



-1000 1000 

Velocity (km/s) 



Fig. 4. Same as Fig.[2]for the remaining 3 galaxies. The final 3 
plots are CO(4-3) spectra from some sources already shown in 
CO(2-l). 



parable to what was found in our lower-z sample. Their average 
is AVfwhm= 370 kms -1 , compared to 348kms~ 1 in paper II. 
Some galaxies have somewhat different widths in their two CO 
lines detected, but this is likely due to the noise intrinsic to the 
data. 



4.2. CO luminosity and H 2 mass 

Since we have not observed the fundamental CO(1-0) line, 
which is a direct measure of the total H2 mass, but mostly CO(2- 
1), it is interesting to compute L'co> the special unit CO luminos- 
ity, through integrating the CO intensity over the velocity profile. 
This luminosity, expressed in units of K km s~'pc 2 , will give the 
same value irrespective of J, if the CO lines are saturated and 
have the same brightness temperature. 
This CO luminosity is given by 



Leo - 23.5/co^b 



(1+z) 3 



K km s pc 



where Iqo is the intensity in K km s , Q# is the area of the main 
beam in square arcseconds and Di is the luminosity distance in 
Mpc. As mentioned above, all sources can be considered to be 
unresolved, as our beam is typically 130kpc in size. We here as- 
sume a ratio of 1 between the CO(2-l) and CO(1-0), luminosi- 
ties (or brightness temperatures) as expected for a warm optically 
thick, and thermally excited medium. In paper II, we found some 
sources with lower excitation, but in the present sample, the exci- 
tation appears higher, as discussed in section l4~3l In any case, our 



s 


Source 


RA(2000) 


DEC(2000) 


z 


SI 


"J005009.8 1-003900.5 


00:50:09.8 


-00:39:01 


0.727 


S2 


FF J0139+0115 


01:39:27.4 


+01:15:52 


0.612 


S3 


IRAS Z02317-0152 


02:34:21.8 


-01:39:01 


0.645 


S4 


"J014612.79+005 112.4 


01:46:12.8 


00:51:12 


0.622 


S5 


IRASZ02433+0110 


02:45:55.5 


+01:23:29 


0.798 


S8 


fl J080430.99+360718.1 


08:04:31.0 


+36:07:18 


0.656 


S9 


"J081929.48+522345.2 


08:19:29.5 


+52:23:45 


0.624 


S10 


"J084846.34+022034.1 


08:48:46.3 


+02:20:34 


0.627 


Sll 


fl J091501.71+241812.1 


09:15:01.7 


+24:18:12 


0.840 


S12 


"J092527. 80+001837.1 


09:25:27.8 


+00:18:37 


0.812 


S13 


"J093741. 65+385752.6 


09:37:41.6 


+38:57:53 


0.621 


S14 


IRAS F10398+3247 


10:42:40.8 


+32:31:31 


0.633 


S16 


SBS 1148+549 


11:51:20.4 


+54:37:33 


0.975 


S18 


3C 280 


12:56:57.1 


+47:20:20 


0.996 


S20 


3C 289 


13:45:26.4 


+49:46:33 


0.967 


S21 


"J144642.28+01 1303.0 


14:46:42.3 


+01:13:03 


0.725 


S22 


"J144953.69+444150.3 


14:49:53.7 


+44:41:50 


0.670 


S23 


IRASF14537+1950 


14:56:04.4 


+ 19:38:46 


0.640 


S25 


"J153244.01+324246.6 


15:32:44.0 


+32:42:47 


0.926 


S26 


"J160222.38+164353.7 


16:02:22.4 


+ 16:43:54 


0.672 


S27 


fl J160825.24+543809.8 


16:08:25.2 


+54:38:10 


0.907 


S28 


CFN1 078 


16:12:37.0 


+54:28:51 


0.910 


S29 


"J161422.14+323403.6 


16:14:22.1 


+32:34:04 


0.710 


S30 


"J172507.39+370932.1 


17:25:07.4 


+37:09:32 


0.689 


S32 


IRAS 17548+6706 


17:54:49.6 


+67:05:57 


0.919 


S34 


IRAS Z21293-0154 


21:31:53.5 


-01:41:43 


0.730 


S35 


PG 0044+030 


00:47:05.9 


+03:19:55 


0.623 


S36 


''EGS70-158 


14:21:37.5 


+53:17:18 


0.930 


S37 


*EGS70-124 


14:22:39.2 


+53:24:21 


0.850 


S38 


"J142301. 15+532419.2 


14:23:01.1 


+53:24:19 


0.660 


S39 


"J142348.80+533006.3 


14:23:48.8 


+53:30:06 


0.770 


SAD 


r_A jo / u-uu / 


It.it.jU. 1 




U.7UU 


S41 


6 EGS70-058 


14:25:02.8 


+53:31:25 


0.690 


S42 


"J143047.32+602304.4 


14:30:47.3 


+60:23:04 


0.607 


S43 


"J143151. 84+324327.9 


14:31:51.8 


+32:43:28 


0.660 


S44 


"J143820.76+340234.3 


14:38:20.7 


+34:02:34 


0.670 


S45 


"J172228.04+601526.0 


17:22:28.0 


+60:15:26 


0.742 


S46 


"J171302.41+593611.1 


17:13:02.4 


+59:36:11 


0.668 


S47 


"J17 13 12.08+600840.4 


17:13:12.0 


+60:08:40 


0.759 



["] SDSS source; [*] [SWR2008]. Redshifts are taken from the NED 
data base. 



hypothesis of a CO(2-1)/CO(1-0) ratio equal to unity can only 
underestimate the CO(1-0) luminosity, and thus the molecular 
masses. Under this assumption, we compute H2 masses using 
M Hj = a L'co, with a = 0.8 M (K kms -1 pc 2 ) -1 , the appro- 
priate factor for ULIRGs. The molecular gas masses are listed 
in Table [3] The choice of a common conversion factor might not 
be realistic, but it is a first approximation before knowing more 
physical quantities about each source, and adapting for each an 
individual factor, based on observation of more CO lines and 
of the gas spatial extent. In the following, CO luminosities are 
often used instead of H2 masses, to remind the reader of this 
uncertainty. 

The average CO luminosity for the 15 galaxies detected is 
L'co = 1.8 10 10 K kms -1 pc 2 , corresponding to an average H2 
mass of 1.5 10 1() M . The star formation efficiency (SFE), also 
listed in Table[3] is defined as L F i R /M(H 2 ) in L /M . The SFR is 
related to the FIR luminosity, where we adopt the relation SFR= 
Lfir /(5.8 1O 9 L ) compiled by Kennicutt (1998); the gas con- 
sumption time-scale can be then derived as t = 5.8/SFE Gyr, 
where SFE is in units of L /M . 
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Table 3. Molecular gas mass, star formation efficiency, dust pa- 
rameters, half-light radius and type 



S M(H 2 ) SFE T rf M d R 1/2 Type 

10" M L Q /M G K 10 s M Q kpc 



SI, S2: Seyfert 1 & 2, Q: QSO, wi, si: weak & strong interaction 
M(H 2 ) and SFE are defined in Sec.l4~8l 



4.3. Molecular gas excitation 

Seven sources have been observed in two CO lines, CO(2-l) and 
CO(4-3), and six have been detected in both, as listed in Table 
[4] One has an upper limit, in the CO(2-l) line. To compare the 
two lines, and derive the excitation, we have computed both the 
total and the peak flux ratio S43/S21 since the CO(4-3) and CO(2- 
1) lines have sometimes different measured linewidths, which 
we attribute to noise. The corresponding ratio between the peak 
brightness temperatures are also displayed in Table [4] to allow 
for easy comparison with the predictions of excitation models. 
The temperature ratio is corrected for the different beam sizes, 
assuming that the sources are unresolved in all our observations. 

The excitation essentially depends on two parameters, the 
H2 volume density, and the kinetic temperature. In ultralumi- 
nous objects, the column densities of the molecular gas are high 
enough (NXFb.) > 10 24 crrT 2 ) that the CO lines are always op- 
tically thick in the low-J lines. This result was reached through 
mapping the gas content of local ULIRGs, where the gas is very 



concentrated (e.g. Solomon etal. 1997). Depending on the ve- 
locity width of the lines, the CO column density per unit veloc- 
ity width (km/s) is larger than 10 17 cirT 2 . To constrain the ki- 
netic temperature of the gas, we have computed the dust temper- 
ature deduced from the far-infrared fluxes (see Table [3), assum- 
ing k v oc V s , where k v is the mass opacity of the dust at frequency 
v, and /3 — 1.5. The average dust temperature for our sample is 
54 ± 5 K. This is higher than what is obtained for local starburst 
galaxies, which have dust temperatures * 40 K (e.g. Sanders 
& Mirabel 1996, Elbaz etal. 2010). Yang etal. (2007) have ob- 
served seven of our objects at 350/im and derive more precise 
temperatures, that are very similar to the values computed above. 
We note a clear increase of the dust temperatures with respect to 
our 0.2 < z < 0.6 sample in paper II. This could be caused by se- 
lection effects, and also due to the frequency used to measure the 
temperature. The 60 and lOO^um bands correspond to rest-frame 
45 and 76/im at z~ 0.32, and rest-frame 35 and 58/vm at z~ 0.72, 
the median redshift of the two samples. 

As in paper II, we assume that the gas is predominantly 
heated by collisions with the dust grains. Indeed, there could 
be small regions that are heated by UV photons of young mas- 
sive stars, or by shocks, but our beam encompasses kpc scales, 
and these would be spatially diluted. The gas kinetic temperature 
should be at most equal to the dust temperature, this is the con- 
straint used in our LVG modeling. We compare the CO-based 
molecular masses and dust masses in section POl 

Using the Radex code (van der Tak etal. 2007), we have 
computed the predicted main beam temperature ratio between 
the CO(4-3) and CO(2-l) lines, for several kinetic temperatures 
and as a function of H2 densities and CO column densities. 
Figure [5] shows these predictions for Tt=60K and 30 K. The 
black contours delineate the range of observed values. In Table|4] 
we list the derived values for the n(H2) densities, for two values 
of the kinetic temperatures (60 and 30K), and for a fixed column 
density per velocity width. 

We adopted a column density of N(CO)/AV of 7xl0 16 
cm" 2 /(kms _1 ), which corresponds to what is expected if the 
high molecular gas masses derived in Table|2]are highly concen- 
trated in the few central kpc. For M(H 2 ) = 3 10 10 M , a typical 
CO abundance of CO/H2 = 10 4 , and a linewidth of 300 kms"\ 
this column density corresponds to a homogenous disk of 3 kpc 
in size. This average column density would be a lower limit, if 
the gas inside 3kpc is clumpy. 

The excitation of the CO gas is in general higher than in our 
lower-z sample. In 4 out of the 7 sources where we have ex- 
citation constraints through measurements of two lines, the CO 
emission ladder should be populated well above the J=4 line. In 
the 3 remaining sources, the excitation is typical of "normal" or 
weakly interacting galaxies, like the Milky Way or the Antennae 
(e.g. Weiss et al. 2007). 

We note that these conclusions on the excitation are prelimi- 
nary. Observations of several lines for more sources are required, 
as there are large variations from source to source. In addition 
the spatial extent of the emission needs to be obtained through 
interferometric measurements. 

4.4. CO luminosity and redshift 

As in paper II, we can now investigate the evolution of the 
molecular gas content of galaxies with redshift. The CO lumi- 
nosity of our sample objects, selected to be bright with high 
masses, is now compared to data reported in the literature at dif- 
ferent redshifts, also selected with large masses, in Figure|6] Our 
full sample (filled black circles, paper II and this work) now fills 
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Log(N c0 /Av) [cm 2 /kms l ] Log(N c0 /Av) [cm 2 /kms l ] 

Fig. 5. Peak Tb ratio between the CO(4-3) and CO(2-l) lines versus H2 density, and the CO column density per unit velocity width 
(Nco/AV) for two values of the kinetic temperature: T* = 60K, the dust temperature (left), and T* = 30K (right). The black contours 
are underlining the values obtained in the data. The predictions come from the LVG models that are part of the Radex code (van der 
Tak2007). 



in the CO redshift desert, between z-0.2 and 1. The average CO 
luminosity over our 33 detected galaxies is L'co = 1.9xl0 1() K 
km s _1 pc 2 , while the average over the local star forming galax- 
ies (Gao& Solomon 2004) is 0.36xl0 10 K km s" 1 pc 2 ,afactor5 
lower. It has been already established that galaxies have more gas 
at high redshift, z ~ 1-2 (Tacconi etal. 2010, Daddi etal. 2010), 
and now it is possible to determine that the rise is gradual, and 
starts already at z=0.2-0.3. This trend will be discussed also for 
the gas fraction (M ?flJ /(M ?flJ +M„) in Section H~6l 

Obviously, such a trend is only indicative and should be con- 
firmed by the study of a large number of galaxies, in the main 
sequence of star formation, and sampling a larger part of the lu- 
minosity function. Here we focus on the most luminous objects, 
and the evolution observed is only that of most actively star- 
forming galaxies. We expect that this evolution in L'co probably 
underestimates the evolution in molecular mass: at higher z, the 
CO luminosity has been measured in higher-J lines, and because 
of likely subthermal excitation, those lines must have lower lu- 
minosity than the CO(1-0) line. 

Theoretically, it is expected that galaxies are richer in molec- 
ular gas at high redshift. Indeed, galaxies at a given mass should 
be smaller (Somerville et al. 2008) which is confirmed by obser- 
vations (e.g. Nagy etal. 2011). As a consequence, the gas com- 
ponent was more compact and thus denser, and the higher gas 
pressure transforms the HI into a molecular gas phase, increasing 
the H2/H1 ratio. Since the abundance of HI is not predicted to 
vary significantly since z=1.5, the H2/H1 ratio reflects the evo- 
lution of the molecular gas content, traced by CO luminosity. 
Obreschkow etal. (2009) and Obreschkow & Rawlings (2009) 
followed the H2/H 1 ratio statistically over 30 million simulated 
galaxies, and predicted its cosmic decline as Qh 2 /^hi 00 (1+z) 1 ' 6 . 
To guide the eye, we have reproduced this behaviour as a red line 
in Figure [6] The amplitude of variation of the L'co envelope ap- 
pears to approximately follow the line. 




log (z) 

Fig. 6. Measured CO luminosities, corrected for amplification 
when known, but not for gas excitation, as a function of redshift. 
We compare our points (filled black circles, and arrows as up- 
per limits) with a compilation of high-z molecular gas surveys, 
and local ones: green triangles are from Gao & Solomon (2004), 
blue squares from Solomon etal. (1997), open circles from 
Chung etal. (2009), blue diamonds from Geach etal. (2009), 
black crosses from Iono etal. (2009), red stars, from Greve 
etal. (2005), green filled circles from Daddi etal. (2010), blue 
asterisks from Genzel etal. (2010), and blue filled circles from 
Solomon & vanden Bout (2005). For illustration purposes only, 
the red curve is the power law in (l+z) L6 for Qh 2 /£2hi proposed 
by Obreschkow & Rawlings (2009). 
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4.5. Correlation between FIR and CO luminosities 



flux Sioo, we can estimate the dust mass as 



It is now well-known that a good correlation exists between FIR 
and CO luminosities (e.g. Young & Scoville 1991). It holds for 
star-forming galaxies, either in the main sequence or in the star- 
bursting phase, and even for quasars, implying that their FIR 
luminosity is dominated by star formation (Iono etal. 2009, 
Xia etal. 2012). The correlation is slightly non-linear, and Xia 
etal. (2012) found a power-law of slope 1.4 as their best-fit, 
which is confirmed in Figure [7] this means that ultra-luminous 
sources are forming stars at higher efficiency, and with a deple- 
tion time as low as ~ 10 Myr. 

In Figure [7] we draw three lines corresponding to three 
orders of magnitude of the ratio Lfir/M(H2), if the CO-to- 
H2 conversion factor adapted for ULIRG is used. Within this 
hypothesis, all galaxies in our sample are located above the 
curve Lfir/M(H2)=100 L /M q (corresponding to a consump- 
tion time-scale of r =58 Myr). However, all detected galaxies 
at high z in the literature are also in this regime, including the 
main-sequence star forming galaxies at high redshift studied by 
Genzel etal. (2010) or Daddi etal. (2010). In those objects, gas 
excitation and spatial extent argue in favor of a more "normal" 
conversion factor, as for the Milky Way. If we adopt a MW-like 
conversion factor, the H2 mass then moves to the right by a fac- 
tor 5.75 (or 0.76 in log), and these objects will fall in the range of 
Lfir/M(H 2 ) between 10 and 100 L o /M with a depletion time- 
scale approaching the Gyr. In such conditions, star formation 
could be sustained continuously, given the frequent cold gas ac- 
cretion from the intergalactic medium. 
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Fig. 7. Correlation between FIR and CO luminosities, for our 
sample (filled black circles, and arrows for upper limits) and 
other galaxies from the literature (same symbols as in Fig[6j). The 
3 lines are for L FIR /M(H 2 )=10, 100 and 1000 L /M from bot- 
tom to top, assuming a conversion factor a = 0.8 M (K kms -1 
pc 2 )" 1 . The three lines correspond to gas depletion time-scales 
of 580 (bottom), 58 (middle) and 5.8 Myr (top). 



In order to better constrain the correct conversion factor to 
apply to each object, we derive the dust mass from the far- 
infrared measurements, and compare it with the gas mass. Given 
the already derived dust temperature T^ and the observed 100 fim 



M(7 



t <l+z)K„B„(T d ) 

5(1 + z)-< 4+ « S l00m {exp( 144(1 + z)/T d ) - 1) 



M„ 



where S vo is the observed FIR flux measured in Jy, £> ] 2 1pc is the 
luminosity distance in Mpc, B vr is the Planck function at the 
rest frequency vr = vo(l + z), and we use a mass opacity co- 
efficient of 25 cm 2 g -1 at rest-frame 100 fim (Hildebrand 1983, 
Dunne etal. 2000, Draine 2003), with a frequency dependence 
of 0=1.5. Estimated dust masses are displayed in Table|3] If we 
adopt the low conversion factor of a = 0.8 M (K kms -1 pc 2 ) -1 , 
the average gas-to-dust mass ratio is 206 for the detected galax- 
ies. The gas-to-dust mass ratio would increase up to 1200 if the 
standard (MW) conversion factor is used. Since the Milky Way 
gas-to-dust ratio is 150, which is also a value typically found in 
nearby galaxies in the SINGS sample (Draine et al. 2007), while 
values larger than 1000 are found only in elliptical galaxies, or in 
low-metallicity dwarfs (Wiklind etal. 1995, Leroy etal. 2011), 
we conclude that the ULIRG conversion factor is not far from 
being adequate for our sample. 



Table 4. CO gas excitation 
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Tb43/Tb21 


n(H 2 ) 


n(H 2 ) 




[total] 


[peak] 
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[cm" 3 ] 


[cm" 3 ] 










60K 


30K 
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2.0 
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1.7 ±0.3 
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n(H 2 ), for T t =60K and 30K, and N co /AV=7xl0 16 cnr 2 /( kms" 1 ) 



4.6. Stellar mass and gas fraction 

The redshift evolution of the CO luminosity envelope found in 
section 14.41 suggests that the declining gas content of galaxies 
could be an important driver of the decline in the cosmic star 
formation density since z=l. To better quantify this, we now 
measure the gas fraction of all systems at different redshift, and 
estimate their stellar masses. 

To compute stellar masses from observed optical and near 
infrared (NIR) magnitudes, standard relations exist as a func- 
tion of colors, derived from stellar populations models, see e.g. 
Bell et al. (2003). We have used them for the local z=0 galaxies. 
For higher redshift objects, K-corrections need to be applied, 
and we used the analytical approximations from Chilingarian 
et al. (2010) in paper II. However, for z>0.5, these are no longer 
valid, and for the present sample we have estimated the stellar 
mass from SED-fitting of the optical and near-infrared lumi- 
nosities taken from public catalogs, mainly the SDSS, 2MASS 
and Spitzer/IRAC fluxes. Most galaxies have only SDSS-DR8 
fluxes, but 10 have in addition 2MASS photometry, and an- 
other 10 have Spitzer-IRAC fluxes. The broadband photome- 
try were fitted using the code FAST (Fitting and Assessment 
of Synthetic Templates) described in Kriek etal. (2009). We 
selected the library of stellar population synthesis models of 



8 



F. Combes et al.: CO in ULIRGs at 0.2 < z < 1.0 



CO 



O 




logSFR=2.24 
Av= 1.6 
logtau= 8.5 

J I I I L_ 



logsSFR=-8.81 
logM*= 11.05 
logAge= 8.6 

J I I I I I ! 



3.5 



4.5 



log(A obs (A)) 



Fig. 8. Spectral energy distribution (SED) of S43 (shown here 
as a representative example), based on optical/near-IR broad- 
band photometry. The red points with error bars represent the 
observations of SDSS (ugriz bands), and Spitzer/IRAC 3.6 and 
4.5/mi observations. The black curve is the best-fit stellar synthe- 
sis model overplotted, obtained with the FAST code (see text). 
An exponential star-formation history, with time-scale tau is as- 
sumed, with a Salpeter IMF. Values for tau and population age 
are given in yr, as well as stellar mass (M ), visual extinction 
(mag) star formation rate (M /yr), or specific SFR (yr -1 ). All 
values are in log, except A v . 
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Fig. 9. The star formation rate (SFR) obtained from the far in- 
frared luminosity, versus stellar mass of galaxies in our sam- 
ple (filled black circles for 0.2 < z < 0.6, and filled blue cir- 
cles for 0.6 < z < 1.0), compared to the sample of Da Cunha 
etal. (2010, filled red squares), Fiolet etal. (2009, filled ma- 
genta circles), Gao & Solomon (2004, filled green triangles), 
Solomon etal. (1997, filled blue squares), Chung etal. (2009, 
open circles), Daddi et al. (2010, filled green circles), and Genzel 
et al. (2010, blue asterisks). The black dashed line represents the 
star forming galaxy main sequence at z=0, the red filled line 
the main sequence at z=l, and the green filled line the main 
sequence at z=2. All three lines have an adopted slope of 0.6 
(Karim etal. (2011). 



Bruzual & Chariot (2003), and adopted an exponential star for- 
mation rate SFR oc exp(-t/tau), with a Salpeter initial mass func- 
tion (IMF, Salpeter 1955). We used the extinction law by Calzetti 
et al. (2000), and metallicity was assumed to be solar. As for the 
IRAC fluxes, only the first two channels (3.6 and 4.5 /mi) are 
relevant for the stellar component at the redshifts under consid- 
eration, and were used for the fit. The visual extinction A,, was 
allowed to vary between and 9, but the best-fits were obtained 
for values always lower than 3. The model fitting gives also an 
estimation of the star formation rates, which is frequently in 
agreement with that derived from the far-infrared luminosities, 
except in a few cases. The time scale for star formation is also 
derived, between the two limits that we set for the model fitting 
10 Myr < tau < 300 Myr. We imposed the latter upper limit, to 
be consistent with the starburst nature of the objects. The param- 
eters of the fit were selected such that the inverse of the specific 
star formation rate is always larger than tau. An example fit for 
the source S43 is presented in Figure[8] There is a certain degen- 
eracy between the physical parameters, since the tau values are 
not tightly constrained, however the stellar mass is more robust, 
depending on the SED fit and magnitudes observed. The derived 
masses are listed in Table [2] 

Stellar masses lie between 4xl0 9 and4xl0 12 M , with a me- 
dian value of 1.1x10 M . Figure [9] displays the star formation 
rate, derived from the infrared luminosity, versus stellar mass, in 
comparison with other samples considered before, and addiing 
the ULIRG samples of Da Cunha et al (2010) and Fiolet et al 
(2009). It is interesting to locate on the graph the position of the 
main sequence of star forming galaxies, as defined by Noeske 



etal. (2007) and Daddi et al. (2007). We adopt the power-laws of 
slope 0.6, and redshift evolution, as SFRoc M° 6 (l+z) 3 5 as found 
by Karim et al. (201 1). Our galaxy points sample the region sig- 
nificantly above the main sequence, with however some scatter, 
so that some galaxies reach the center of the main sequence. The 
star-forming galaxies at z=0, whose FIR luminosities are not all 
ULIRG (green triangles) are also somewhat above the z=0 main 
sequence. 

The gas fractions F,, as =MH2/(M(H2)+M») derived from our 
stellar masses show large variations (see Figure ITOt. To obtain 
the gas fractions, we have converted the CO luminosity to a gas 
mass using the ULIRG factor (o , =0.8). For the high-z samples 
(blue asterisks and green dots), where the standard MW conver- 
sion has been selected by the authors, we have in addition ex- 
tended the points by a dotted line joining the two extreme values 
of gas fraction, obtained with <?=0.8 and 4.6 (the MW value). 
The gas fraction is well correlated with CO luminosity. The av- 
erage gas fraction in our galaxies is 15% for the 0.2<z<0.4 sam- 
ple, and 24% for the 0.6<z<1.0 sample. For the low-z samples 
by Gao & Solomon (2004), Solomon etal. (1997) and Chung 
et al. (2009) the averages are 5%, 6% and 7%, respectively, while 
they are 23% and 10% for the samples of Daddi etal. (2010) 
at z=l .5 and Genzel etal. (2010), at z=l-2. Assuming a stan- 
dard MW conversion factor, the two last values become 63% 
and 35%. The increase of gas fraction with z is clearly apparent, 
as shown in the Figure [TT] To guide the eye, we have indicated 
in this figure the z-evolution of the cosmic star formation his- 
tory, as compiled by Hopkins & Beacom (2006), from different 
works in the literature, and complemented at very high redshift 
by the gamma-ray burst (GRB) data of Kistler et al. (2009) and 
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the optical data (Lyman-Break Galaxies, LBG) from Bouwens 
etal. (2008). In this approximate comparison, it is interesting to 
note that the rise of the SFR has only slightly higher amplitude 
than the rise of gas fraction, supporting the hypothesis of the 
large role of the gas in this evolution. 
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Fig. 10. Gas fraction=MH 2 /(M(H2)+M,,), versus CO luminosity, 
assuming the same CO-to-H2 conversion factor a = 0.8 M (K 
kms~' pc 2 )~'. The points of the high-z samples (blue asterisks 
and green dots), have been continued by a dotted line joining the 
two extreme values of gas fraction, obtained with ff=0.8 and 4.6. 
All symbols are as defined in Fig [6] 
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Fig. 11. The gas fraction as in Figure [TO] but versus redshift z. 
The red curve is a schematic line summarizing the evolution of 
cosmic star formation density, from the compilation by Hopkins 
& Beacom (2006), complemented with the GRB data by Kistler 
et al. (2009) and the optical data from Bouwens et al. (2008). The 
red curve is logarithmic and only indicative of relative variations 
of the star formation rate per cubic Mpc as a function of redshift, 
and can be translated vertically. 



4. 7. Activity of the galaxies 

In order to determine what is triggering and regulating the star- 
burst activity in our sample galaxies, it is interesting to try to 
determine which galaxies are morphologically perturbed or in- 
teracting, and to infer the activity status of their nucleus. The last 
column of Table [3] indicates the classification of nuclear AGN 
activity and/or possible minor or major interactions. The latter 
have been derived from the available images, either from HST, 
or presented by Stanford et al. (2000), or SDSS images. Some of 
the HST images of galaxies detected in CO are shown in Figure 
IT2l the remaining in the Appendix. The strong or weak interac- 
tions were determined from the relative fraction of the light in 
the tidal tails or perturbed features. No interaction means that 
the galaxy image looks unperturbed. 

At large infrared luminosities, the fraction of Seyfert is about 
50% in local ULIRGs, and their fraction is increasing steeply 
with L FIR (Veilleux etal. 1999). The tight link between CO and 
infrared luminosities however shows that the dust heating is 
dominated by the starburtst in these objects (Iono etal. 2009). 
Among our 39 objects, 13 are known as AGN, and 7 of these 
are detected in CO, implying a detection rate of more than 50%, 
higher than that of the whole sample. The percentage of weak 
and strong interactions in the whole sample are 23% and 28%, 
respectively, while in the detected objects sample they are both 
33%. There is thus a slight correlation between CO detection and 
interactions, however not at high significance. For more stringent 
conclusions, higher spatial resolution images need to be obtained 
for the half of the sample without HST data. 



4.8. Star formation efficiency 

Like in paper II, we define the star formation efficiency as SFE 
= Lfir/M(H2), assuming a constant CO-to-H2 conversion fac- 
tor. The average SFE in our 0.6 < z < 1 sample is 595 L /M o , 
quite comparable with that of the 0.2<z<0.6 sample. We plot 
SFE versus Lfir in Figure Qj] and versus redshift in Figure 1X41 

Our intermediate-z sample shows among the largest efficien- 
cies in star formation. The most extreme objects, with SFE > 
1000 L /M are mostly interacting or perturbed. While there is 
a relatively good correlation between SFE and Lfir, it is not the 
case for L'co- The latter shows only two populations of objects, 
with the local starbursts at low efficiency and gas content, sepa- 
rated from the higher-z samples. The SFE correlates much better 
with the dust temperature; this was also found in paper II, and 
we do not reproduce the figure here. 

Finally, the evolution of SFE with redshift is shown in Figure 
[T4l There is an obvious rise of the envelope between z=0.2 and 
1, precisely in the range of our sample. The amplitude of the rise 
appears however lower than that of the cosmic star formation 
history, as compiled by Hopkins & Beacom (2006), and shown 
here as the red line. Figure [TT] has shown that the redshift vari- 
ations of the gas fraction have a larger amplitude. The star for- 
mation evolution is certainly due to a combination of factors, 
essentially the gas fraction and the efficiency. Both vary with 
redshift, but it can be concluded that the evolution of gas frac- 
tion appears as the most determinant parameter, followed by the 
star formation efficiency. 

To search for a possible correlation between the SFE and 
the compactness of the starburst, we have computed the half- 
light radii of galaxies in our sample, using the GALFIT software 
(Peng et al. 2002). Red images in the I-band were used at in- 
termediate z, with the highest possible angular resolution: HST 
images when available, and some with SDSS images. A single- 
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Fig. 12. Optical or NIR images from HST. Most images are ACS-F814W, except S35 (WFPC2-702W) and S47 (WFC3-160W). All 
sources are detected in CO. Each panel is 5"x5" in size, and is centered on the galaxy coordinates given in TableQ] The brightness 
scale is logarithmic (except for S47). North is up and East to the left in all panels. S14, S29, S34 and S47 were classified strongly 
interacting, while SI and S35 are weakly interacting. 
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Fig. 13. Star formation efficiency SFE=Lfir/M(H2), versus far- 
infrared luminosity Lfir, assuming a single CO-to-H2 conver- 
sion factor of a = 0.8 M (K kms~' pc 2 )~'. All symbols are as 
defined in Fig [6] 
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Fig. 14. Same as FigureQj] but versus redshift. The red dash line 
is the same as in Figure [TT] The SFE shows much less evolution 
with z than gas fraction. 



component Sersic model was fit to all images, giving an effec- 
tive radius, which is displayed in Table [3] We compare to the effective radii in the B-band for some local galaxies, obtained in 
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HyperLeda (Paturel et al. 2003fl The SFE is plotted versus the 
half-light radii of galaxies R1/2, when available, in Figure [15] 
There is only a slight trend for SFE to be anticorrelated to R1/2, 
or correlated to the compactness of the blue light distribution, 
which reflects the star formation. The compactness of the star- 
burst might be better traced by the molecular gas distribution, 
but this requires high resolution mapping of the high-z galaxies. 
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Fig. 16. Evolution with redshift of averaged quantities, SFE at 
left, and gas to stellar mass ratio at right. The average of detected 
points only is plotted in green, and with the 3<x upper limits in 
black (for high-z samples only). The error bars are the statistical 
ones following the square root of the number of points averaged. 
The red line is the same as in Figures[TT1and[T4l 



CO(2-l) and CO(4-3) transitions to derive the gas excitation. 
The latter varies from source to source, but we find that in the 
majority the gas is significantly excited up to J=4, suggesting a 
high H2 volume density and/or temperature. The comparison be- 
tween the CO-derived gas mass and the dust mass derived from 
the far infrared fluxes supports the choice of the ULIRG conver- 
sion factor for our sources. This does not exclude that some of 
the sources have a low excitation, and their total gas mass has 
been under-estimated. Spatial information on the CO emission 
is needed to better constrain this issue, and we will report the 
results of PdB interferometer imaging in a subsequent paper. 



5. Discussion and conclusions 

We have presented the completion of our survey of the gas con- 
tent in galaxies at intermediate redshift, through CO observa- 
tions in a sample of 39 ULIRGs with redshifts between 0.6 and 
1.0. Together with our earlier work (paper II), this eliminates 
the CO desert in this highly relevant epoch (0.2<z<l), in which 
the star formation rate in the universe decreased by a factor 10. 
The CO detection rate between z=0.6 and 1 is 38%, significantly 
lower than between z=0.2 and 0.6 (60%, cf paper II). The CO lu- 
minosity, and therefore the derived H2 mass are increasing with 
redshift by about a factor 4 up to z= 1 . We have estimated the stel- 
lar mass of all sample objects through SED fitting of the optical 
and near-infrared fluxes, and derived gas fractions. The evolu- 
tion of the gas fraction with redshift is very pronounced, and is 
reminiscent of the star formation rate evolution, suggesting that 
the gas fraction has a large role in determining the star formation 
history. 

We have also checked the evolution of star formation effi- 
ciency with redshift, and found less of a trend compared to the 
gas fraction. We conclude that both parameters have a role in 
setting the SFH, but the gas fraction is more important. This can 
be seen in the envelope traced by the most extreme objects, but 
also in the averages made over all detected objects, with or with- 
out taking into account the upper limits (Figure [16). We have 
checked whether the SFE is related to the compactness of the 
starburst, but found only a slight trend of SFE to decrease with 
half-light radius (Figure \15[. 

Since the gas fraction depends on the CO-to-H2 conversion 
ratio adopted, we observed a number of galaxies in both the 

1 http://leda.univ-lyonl.fr 
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Fig. 17. The gas fraction as in Figure[10j but versus stellar mass. 
The red line represents the predictions of cosmological simu- 
lations at z=l, and the black dashed line at z=0 from Dave 
etal. (2011). 

Several physical parameters may intervene to account for 
the evolution of galaxies since z=l. First the gas fraction is a 
key parameter, and we have shown that indeed, its evolution in 
this redshift range is significant. There is also the star forma- 
tion efficiency, which could be higher, being triggered either by 
galaxy interactions and mergers, which have been more frequent 
in the past, but also through cold gas accretion, which is also 
thought to be efficient at these epochs (Keres et al. 2005). The 
latter mechanism can also produce high gas fraction, resulting 
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in instability-driven turbulence,perturbed disks, and clumpy gas 
distributions. This could mimick galaxy interactions in the ob- 
servable morphology. Between z= 1 and 0, hydrodynamical and 
semi-analytic studies predict a quenching of star formation, be- 
ginning at high stellar mass, and progressively involving lower 
masses (e.g. Gabor etal. 2010, Dave etal. 2011). It is interest- 
ing to compare observations of the gas fraction versus mass, to 
better constrain the models, since they have not yet reached co- 
herence with observations. Figure[17]displays this relation, with 
the model lines superposed for z=0 and 1. They represent the 
best-fit models, corresponding to the no-wind simulation from 
Dave etal. (2011). The consideration of stellar winds will pro- 
vide a quenching mechanism, but leads to an under-prediction of 
the gas fraction in small galaxies. Other quenching mechanisms 
are then necessary for massive galaxies, such as the influence of 
AGN feedback. 
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Fig. 18. The star formation efficiency SFE (which varies as the 
inverse of the depletion time-scale for the molecular gas) versus 
the specific star formation rate sSFR, in yr . The red line is the 
fit from the COLD GASS sample of star forming galaxies at z=0, 
from Saintonge et al. (201 lb). 

Another interesting parameter is the specific star formation 
rate, or SFR normalized by the stellar mass, which gives the 
time-scale of formation of all the stellar mass in a galaxy, at the 
given present SFR. We observe a good correlation with the star 
formation efficiency, in Figure[l8] Note however that the quanti- 
ties on the vertical and horizontal axes are not independent, since 
they both contain SFR. This correlation was studied in particular 
by Saintonge et al. (201 lb), interpreting the SFE as the depletion 
time scale for the molecular gas. With respect to our definition 
of SFE, the depletion time-scale ta ep in yr is such that log(t^ e; ,) = 
9.76 -log(SFE). The diagram helps to identify clearly the main 
sequence of star formation, where galaxies have a continuous 
SFR, and the depletion time-scale should be comparable to the 
stellar-mass forming time-scale. We have superposed the best 
fit line for the sample COLD GASS of 222 normal star forming 
galaxies at z=0. Our galaxies are generally below the line, which 
is due to their larger gas content. 
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Appendix A: Available images for more sources 

We present in this Appendix the remaining HST images of the 
sources, which give insight in their morphology. 
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Fig. A.l. Optical red (F814W) HST-ACS images of six more sources. Only S16 is detected in the CO line. Each panel is 5"x5" in 
size, and is centered on the galaxy coordinates of Table Q] The brightness scale is logarithmic. North is up and East to the left in all 
panels. 




Fig. A.2. Same as in Fig|ATj except for S18, which is a NICMOS-F160W image. S22 and S30 are detected in CO. 
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